Spin excitations in the antiferromagnet NaNiC>2 
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In NaNi02 , Ni 3+ ions form a quasi two dimensional triangular lattice of S = 1/2 spins. The 
magnetic order observed below 20K has been described as an A type antiferromagnet with ferro- 
magnetic layers weakly coupled antiferromagnetically. We studied the magnetic excitations with 
the electron spin resonance for frequencies 1-20 cm -1 , in magnetic fields up to 14 T. The bulk of 
the results are interpreted in terms of a phenomenological model involving bi-axial anisotropy for 
the spins: a strong easy-plane term, and a weaker anisotropy within the plane. The direction of 
the easy plane is constrained by the collective Jahn- Teller distortion occurring in this material at 
480 K. 

PACS numbers: 71.27+a, 75.30.Et, 71.70-d, 61.10Nz, 75.40Cx 
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I. INTRODUCTION 



A two-dimensional triangular network of magnetic ions 
interacting via an antifcrromagnetic interaction is a well 
known geometrically frustrated system where unconven- 
tional magnetic properties are expected 0. Usually, a 
long range magnetic order occurs at low enough tempera- 
ture. For instance, in XCI2 with X=Cr, Br, or in ACr02 
with A=Li, Ag, Cu, or CsCuCl3 the magnetic order is 
based on a 120° spin structure on the triangles. All these 
compounds have an easy plane or easy axis anisotropy 
together with Heisenberg type antiferromagnetic inter- 
actions. In other compounds, no magnetic order was de- 
tected so far (NaCr0 2 , KCr0 2 , NaTi0 2 , LiM0 2 ). The 
possibility of the orbital order competing with the spin 
order makes LiNi0 2 particularly interesting. In this com- 
pounds the Ni 3+ ions have a spin 5 = 1/2 and the or- 
bital occupation is doubly degenerate: the e g orbitals, 
I 3z 2 — r 2 > and | x 2 — y 2 >, have the same energy unless 
the oxygen octahedron around the magnetic ion becomes 
elongated or compressed due to the Jahn- Teller effect. 

Surprisingly, no orbital order has been observed in 
LiNi0 2 . The absence of both orbital and magnetic or- 
der in this compound has been the subject of intense 
debate lately, both experimentally 0] and theoretically 
0. In particular, there is still a controversy on what are 
the relevant magnetic interactions within the triangular 
planes and to what extend this compound is magnetically 
frustrated. The comparison with the isomorphic NaNi0 2 
is aimed at elucidating this unconventional behavior. In 
NaNi0 2 a ferro-distortive orbital order (a collective Jahn- 
Teller distortion) is observed below 480 K |4J and a loa 
range antiferromagnetic order appears below 20 K 
This magnetic order was first described as an A type 



antiferromagnet with ferromagnetic planes coupled an- 
tiferromagnetically 5] 6]. The magnetic superlattice has 
indeed been observed recently in neutron diffraction mea- 
surements 0- However, the description of the magnetic 
system with just two magnetic interactions (an antifcr- 
romagnetic Jaf between subsequent NiO planes and a 
ferromagnetic Jf within the planes) fails to describe the 
whole magnetic behavior, in particular the presence of 
three characteristic fields observed in the magnetization 
curve The direction of the spins measured by neu- 
tron diffraction is also unusual: they point toward the 
center of one of the triangles in the oxygen octahedron 
surrounding each Ni ion, at 100° from the Ni plane 0. 
Moreover, a recent inelastic neutron study has shown 
that a gap of ~0.7 meV is present in the spin excitation 
spectrum. 



We have performed Electron Spin Resonance measure- 
ments in magnetic fields up to 14 T on powder samples 
of NaNi0 2 . We adopted and extended the model used 
to describe the magnon spectrum We find that the 
low field behavior of the system is characteristic of an 
"easy-plane" magnet, with a small anisotropy within the 
plane. The spin-flop transition is assigned to this latter 
anisotropy. At high fields all spins are aligned parallel 
and the saturation effects dominate the behavior. Al- 
though the model described here is generally successful, 
we also observed a spin resonance mode that remains 
unexplained. The possible implications for the magnetic 
interactions in LiNi0 2 are discussed. 
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FIG. 1: Shift of the resonance lines position from the param- 
agnetic line at 75 GHz as a function of temperature. Insert: 
typical ESR spectrum as a function of magnetic field recorded 
at 2.4 K.The dotted line corresponds to the paramagnetic line 
at g=2 . 



II. EXPERIMENTS 

The NaNiC>2 powdered sample was obtained following 
the procedure described elsewhere 0. The monoclinic 
C2/m crystal structure was checked by X-ray powder 
diffraction. The system is quasi two-dimensional: The 
Ni-Ni interplane distance is large, 5.568 A. Below the 
Jahn- Teller transition at 480 K the Ni 3+ ions are ar- 
ranged in a slightly distorted triangular network, with 
one short length and two longer ones (at room tempera- 
ture, the Ni-Ni distances are 2.84 A and 3.01 A respec- 
tively compared to 2.96 A at 565 K) Q. Also, the oxy- 
gen octahedra surrounding the nickel ions are elongated. 
Such a distortion favors the | 3z 2 — r 2 > orbitals. Indeed, 
the ESR spectra at 200K have the anisotropic shape typ- 
ical of the | 3z 2 — r 2 > orbitals, with gW = 2.03 and 
g 1 - = 2.28 for fields parallel and perpendicular to z, re- 
spectively m. 

The magnetic behavior of our NaNiC>2 sample was 
the same as the one reported previously |9j. The mag- 
netic susceptibility exhibits a Curie- Weiss behavior above 
100 K. The Curie- Weiss temperature of Tcw=36 K re- 
flects the predominance of ferromagnetic interactions. 
This is a noticeable difference compared to frustrated 
triangular compounds such as LiCr02- The effective 
magnetic moment (1.85 /is) is in agreement with the 
| 3z 2 — r 2 > orbital configuration for the e g electron in 
the low spin state of Ni 3+ (spin S=l/2, average gyromag- 
netic factor g av = 2.14). The Neel temperature, observed 
as a peak in the susceptibility, is Tn=20 K. The magne- 
tization curve in the ordered phase at 4 K presents two 
kinks at Hco—l-8 T, i?ci= 8 T, previously assigned to 
spin flop transitions, and the magnetization is fully sat- 
urated at H sat = 13 T 0. 

The Electron Spin Resonance measurements (ESR) 



were performed using two different techniques. Classi- 
cal high frequency, high magnetic field ESR was used at 
fixed frequencies using Gunn diode oscillators and mul- 
tipliers at 35, 75, 150 and 225 GHz, combined with a 
9 T magnet at the Budapest University of Technology 
and Economics. Previous measurements at the Grenoble 
High Magnetic Field laboratory have been described in 
Ref. H 

Other measurements were performed by a novel 
method [ToL ITl| . We used the far infrared facilities of the 
National Synchrotron Light Source in Brookhaven Na- 
tional Laboratory, at the U12 IR beamline, including a 
14 T superconducting magnet (Oxford Instruments), and 
a Fourier Transform Infrared Spectrometer (Sciencetech, 
SPS 200). We will present results obtained in the fre- 
quency range 3 cm" 1 - 20 cm" 1 (90 GHz - 600 GHz) at 
4 K. The powdered sample was in a disk-shaped teflon 
container of 5 mm diameter and 1 mm thickness, in the 
center of the magnet. Light pipes were used to guide the 
infrared light from spectrometer to the sample, and the 
transmitted light was detected with a bolometer operat- 
ing at 1.2K temperature. The light propagated parallel to 
the static magnetic field. Spectra were recorded at fixed 
fields. The frequency resolution was selected at 0.2 cm" 1 , 
much less than the width of the resonance lines. The 
upper frequency cut-off was adjusted to about 20 cm" 1 
with a Fluorogold filter of appropriate thickness. The 
lower cut-off was limited by the incident spectrum and 
the spectrometer performance. The signal to noise ratio 
was good down to 3 cm -1 , but there were strong minima 
in the incident light intensity around 3.7 and 4.6 cm" 1 . 
We eliminated the experimental points at these minima 
without loosing the data around them, obtaining reliable 
signal down to 3 cm" 1 . After averaging over several spec- 
tra (typically between two and four) and smoothing, the 
absorption spectrum at a given field was divided by a 
reference spectrum recorded at zero magnetic field, fol- 
lowing the same procedure as described in Ref. [Tlj . 



III. RESULTS 

First, we present the ESR results obtained with the 
classical ESR technique. A typical spectrum (the deriva- 
tive of the absorbtion signal), recorded at 2.4 K at 
75 GHz in the antiferromagnetic phase, has three res- 
onance features (Fig. The lines labeled (A), (B), 
and (C) are shifted from the paramagnetic g = 2 line 
differently and they exhibit different temperature depen- 
dencies. Mode (A) has been already studied in |{| and 
was related to the spin flop field at 1.8 T. Mode (B) re- 
mains close to the paramagnetic line. Mode (C) consists 
of two close resonances which broaden with temperature 
and are hardly resolved above 6 K. Therefore only the 
first, larger intensity, resonance is plotted for the tem- 
perature dependence of mode (C). Below the phase tran- 
sition temperature Tn, this mode behaves as the order 
parameter. 
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The paramagnetic line splitting at 200 K (empty sym- 
bols in Fig. 1) is in agreement with the line splitting 
studied in detail at several frequencies by Chappel et al. 
PJ . There is also an agreement in the temperature depen- 
dence of the splitting. We will discuss the paramagnetic 
state in a separate publication |Tfij . 

The measurements were done on powder samples, and 
therefore the large line broadening seen in Fig. \T\is not 
surprising. For any given crystallite in the sample, the 
line position depends on the relative direction of exter- 
nal field and the crystallographic axes. Relatively sharp 
features appear in the spectrum because in the powder 
sample certain resonance frequencies (e.g. the highest 
and the lowest possible values) acquire large statistical 
weight. This effect is well known from the powder sam- 
ples with significant <?-factor anisotropy. 

It is also important to notice that the presence of three 
lines in the spectrum at finite field does not necessarily 
mean that the microscopic Hamiltonian of the system 
has three modes, for the same reasons as above. On the 
other hand, the number of modes seen at zero external 
field (where all crystallites are equivalent) are directly 
relevant to the Hamiltonian, as we will discuss later. 

The low temperature ESR was further explored by far 
infrared spectroscopic methods. The advantage of detect- 
ing ESR by far-IR spectroscopy is that one can readily 
map the power absorption over the full range of magnetic 
fields and frequencies. The measurement does not rely on 
sweeping the magnetic field; this becomes especially im- 
portant when one tries to discern features at zero field, 
or compare ESR to inelastic neutron scattering done in 
zero field. 

Typical absorbtion spectra at T and 14 T recorded at 
4 K are given in Fig. The sharp peaks in the spectra 
are "instrumental": either due to interference features 
in the synchrotron source, or generated by the multiple 
reflections in the sample, or in the spectrometer itself. 
Nevertheless, the absorption due to spins is clearly vis- 
ible with features around 6 cm" 1 (labeled (1)), 9 cm" 1 
(labeled (2)) and 11-18 cm" 1 (labeled (M)). 

The key to these measurements is taking the ratio of 
two data sets at different fields, eliminating the features 
in the raw spectrum of the incident light, and leaving 
field dependent part of the absorption. We have chosen 
always the same reference: the spectrum at zero field, so 
that the relative transmission J r is related to the sample 
absorption I sp i n (H): 
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I r {H)=I N (H)/I spin (0). 



(1) 



The relative transmission for H=1A T is plotted in the 
bottom part of fig. [5] The data below 3 cm" 1 are not re- 
liable. A smoothing procedure using a FT filter has been 
used to improve the signal to noise ratio. It is responsi- 
ble for the weak oscillations with the 2 cm -1 periodicity. 
The relative transmission data for several magnetic fields 
are plotted in Fig. [31 for the frequency range 4 cm" 1 -20 
cm" 1 . An average over 3 to 4 data has been used. Be- 
low 10 cm -1 two broad peaks ((1) and (2)) lie above 1, 
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FIG. 2: upper part: Typical FT-ESR spectra recorded at 4 
K, T and 14 T. Lower part: Intensity ratio. 
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FIG. 3: FT-ESR Relative transmission I r (H) as a function 
of the electromagnetic wave frequency at different magnetic 
field values H. 



indicating that the absorption / sp i n (0) presents two reso- 
nance modes. From all these relative transmission data, 
we can deduce the frequency dependence of J sp m(0) and 
then I S pi n (H) = I r (H)xI sp i n (0) for each magnetic field 
value. 

All these results can be mapped on a frequency / mag- 
netic field diagram (Fig. 0J. The symbols correspond 
to maxima in the derivative of the ESR absorption mea- 
sured by field sweep at fixed frequencies, the colors are 
derived from the absorption I S p- m (H) measured at fixed 
fields. The main branch (M) is a broad signal, with fre- 
quency approximately proportional to the field, signifi- 
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FIG. 4: Magnetic modes at 4 K in a frequency versus field 
map obtained from the resonance lines positions in conven- 
tional ESR(closed symbols: this work; open symbols: from 
H) and the transmission J S pm in FT-ESR (the color scale on 
the top refers to the intensity of the signal, as shown in Fig. 

The dash black line corresponds to the paramagnetic res- 
onance at g — 2. The blue and red continuous lines represent 
the calculated frequencies oj" and ui ± corresponding to fields 
parallel and perpendicular to the hard axis, respectively, the 
red dotted lines, uj mt , are for fields along the intermediate 
axis (see appendix B). 



cantly above the g = 2 line. Mode (A) and (B) are re- 
lated to the spin- flop transition observed at 1.8 T in the 
magnetization curve @. Mode (A) has been discussed 
previously [f| in terms of an easy axis antiferromagnet 
model p| . The field dependence of this mode is typical 
of crystallites in the powder where the static field is per- 
pendicular to the easy axis. When the magnetic field is 
along the easy axis, a spin-flop transition occurs at 1.8 T. 
Crystallites oriented in this direction contribute to Mode 
(B). 

Modes (1) and (2) are seen mostly in the far-IR mea- 
surements, with supporting evidence from a field sweep at 
225 GHz. These modes have very little field dependence. 
The corresponding zero field gaps are 6.5 and 9.0 cm , 
respectively (0.85 and 1.1 meV). Mode (1) at 0.85 meV 
is close to the dispersionless magnon mode observed at 
0.7 meV in inelastic neutron scattering measurements 0. 
Mode (1) crosses the main branch, and it appears again 
at low frequencies, labelled as mode (C), also seen in 
Fig. ^ Apparently, this mode softens to zero frequency 
around 8 T. 

Further information can be deduced from a quantita- 
tive analysis of the IR absorption data. We will assume 
that the absorbtion I spin (H) can be approximately re- 
lated to the dissipative spin susceptibility x" : 



In 7 spin (H) 



-Airnd 



(2) 



where c is the vacuum speed of light, n the refraction 
index, d is thickness of the sample and uj the light fre- 
quency. The absorption peaks x" (H) are well fitted with 
a Gaussian lineshape (see insert in Fig. Eb). The corre- 
sponding resonance frequency, linewidth and total area 
are plotted as a function of the applied magnetic field 
(FigEt, b ,and c respectively). Note that the line po- 
sitions are unchanged whether one looks at I sp in(H) or 
X" (H) . The line position and line width confirm that the 
resonance associated with the main branch (M) is nearly 
paramagnetic with g « 2.2. Its linewidth is in accordance 
with the linewidth calculated from the g anisotropy at 
200 K. However the anisotropic line shape observed at 
200 K is no longer present here, presumably due to addi- 
tional lifetime broadening. The branches (1) and (2) are 
clearly identified at 6.5 cm -1 and 9 cm -1 with no field 
dependence until they both merge into the main branch 
(M). 

The area under the absorption peak is related to the 
ESR susceptibility through the Kramers-Kronig relation: 



X (H,0) cx f ^ H,w) du) ~ — / x"(H,oj)dw (3) 

J u m J 

where ujq is the peak position. For the main branch (M), 
uj is proportional to the static field H. Then the peak 
area J x"(H,u))du> scales with the static magnetization. 
This is indeed observed in Fig. [3J;. The saturation of 
the total intensity above 9 T is in accordance with the 
saturation of the static magnetization observed in this 
material 0|. For branches (1) and (2), wo is constant. 
Then the peak area is proportional to the static suscep- 
tibility x(H,0). For both modes, it is constant with the 
field, as is expected away from spin-flop field transitions. 
One can also see that mode (2) has a weaker suscepti- 
bility than mode (1). Note that this susceptibility mea- 
sured on a powder sample includes the statistical weight 
associated with the crystallites in the powder which are 
oriented correctly as regards the magnetic field direction 
for each resonance mode. 



IV. DISCUSSION 

Several models have been proposed to describe the 
magnetic properties in this material. The first one as- 
sumes a uniaxial easy axis antiferromagnet, proposed by 
Bongers et al. [fj to account for the spin-flop transi- 
tion at 1.8 T and then used in Ref. Q for the ESR 
magnon branch at 53 GHz (branch (A) in Fig. QJ. This 
model includes a strong ferromagnetic coupling in the Ni 
planes, a weaker antiferromagnetic coupling between the 
Ni planes and an even weaker anisotropy which aligns 
the spins in the easy axis direction. It was shown later 
that this model was incomplete in terms of explaining 
the magnetic properties, and an additional energy scale 
was required [Jj. This conclusion was confirmed by re- 
cent neutron inelastic scattering measurements 8j. The 
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FIG. 5: Gaussian fit of FT-ESR susceptibility at 4 K: line po- 
sition (a), line width (b)and line total area (c) together with 
the sample magnetization (continuous line) as a function of 
applied magnetic field. The dotted line in (a) corresponds to 
the paramagnetic line at g=2.2 and in (b) to a g anisotropy 
linewidth A<7=0.26. Insert: gaussian fit of the ESR suscepti- 
bility. 



magnon spectrum was interpreted in terms of an easy 
plane AF model with the same ferromagnetic and anti- 
ferromagnetic couplings as in the previous model. The 
anisotropy was assigned to the ferromagnetic exchange 
coupling, and its magnitude was an order of magnitude 
larger than the weak anisotropy used in the evaluation of 
the ESR data. 

These models were merged by assuming an anisotropy 



tensor that has an easy axis within an easy plane, with 
a much larger energy needed for turning the spins out of 
the plane than within the plane. This simple approach 
accounts for most of the experimental observations, in- 
cluding a better description of the low field ESR data of 
Chappel et al. . We start with the following anisotropic 
spin Hamiltonian: 



Ho 



Sj+gfisSU (4) 



where H is the magnetic field, ^2 ab is over nearest neigh- 
bor spins within the layers, and J] is over nearest neigh- 
bor spins along the c directions. The exchange interac- 
tions are represented by tensors: J a b describes the ferro- 
magnetic intraplane exchange, and J c is the interplane 
antifcrromagnetic exchange. The intensity and the fre- 
quency of the spin resonance lines for a crystal positioned 
in arbitrary direction relative to the external field can be 
calculated from the Hamiltonian in two steps. First, the 
spin configuration is determined by minimizing the free 
energy; second the ESR frequencies are related to the 
small oscillations around the equilibrium configuration. 
There are two resonance modes for any given direction 
when the magnetic order is described by two sublattices. 
Occasionally (in zero field or field applied along high- 
symmetry directions) the two modes may be degenerate. 

The powder spectrum is obtained by the average of 
the absorption of individual grains over all directions. 
In practice, the calculation of frequency and intensity 
for arbitrary field orientation is rather hard. Instead, 
the resonance frequency is calculated for the static field 
pointing in three principal directions. The six "principal 
frequencies" (three frequencies for each of the two modes) 
play a special role in the powder average, since some of 
them correspond to extremal values and the statistical 
weight factor for these frequencies will be high. 

This process is similar to the determination of the *g* 
tensor from the paramagnetic resonance of a powder sam- 
ple with a g-factor anisotropy. Naturally, the orientation 
of the principal axes of the tensor relative to the crystal- 
lographic directions cannot be determined from measure- 
ments on a powder sample. Nevertheless, the ^-tensor 
components are routinely determined from powder mea- 
surements this way. It is often sufficient to look at the 
highest and lowest frequencies in the measured broad 
ESR line, and one can identify the highest and lowest 
tensor components. 

There is an extensive literature on the anisotropy ef- 
fects in antiferromagnets, starting with the early work of 
Keffer and Kittel nJ on the easy axis problem. The 
three axis case (easy, intermediate, hard) was treated 
by Nagamiya Here we base our analysis on this 

work and on a review by Turov |14| . The magnon spec- 
trum is treated in the linear approximation, equivalent to 
a quasi-classical treatment and we focus on the ground 
state properties (T=0 K). The details of the calculation 
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are given in Appendix A and B. First, the Hamiltonian 
is transformed into: 

Hmf/N = ^m 2 + + + ^ + _ mhj (5) 

where m and 1 are the total magnetization and the anti- 
ferromagnetic order parameter, respectively. Parameter 
A describes an isotropic antifcrromagnetic coupling, a 
and c describe the anisotropy of the total ferromagnetic 
moment. In the ideal antiferromagnetic state m = 0, 
therefore these parameters have no influence on the free 
energy of the system in zero external magnetic field. On 
the other hand, parameters b and d act on the order 
parameter 1 directly. Kittel treated the uniaxial model 
(c — d = 0) at low fields; easy plane and easy axis cor- 
responds to b > and b < 0, respectively. Nagamiya's 
biaxial model is obtained when a = c = 0. Turov dis- 
cusses the uniaxial case for fields up to the saturation 
field [13 ■ 

We set b > so that the zero field equilibrium mag- 
netizations are in the plane perpendicular to z; this will 
be the "easy plane". The anisotropy within the plane, 
represented here by coefficients c and d, is much smaller: 
c,d << b. (Notice that these choices do not follow from 
the symmetry arguments related to the lattice distortion, 
but they were forced by the experimental observations). 
Since c and d have similar effects on the spin resonance 
frequencies, and we can not determine them separately, 
we will assume that d = 0. This leaves us with four pa- 
rameters to determine: A, a, b and c. For c < 0, x will be 
the easy direction, y will be the intermediate direction, 
and z remains the hard direction. 

To facilitate the discussion, it is convenient to intro- 
duce the "effective fields" H E = A/M , H' E = (A + 
a)/Mo, and H' EA — y/Ac/Mo where Mo is the saturation 
magnetization (more effective fields follow in Appendix 
B). For applied field perpendicular to the hard axis the 
spin system saturates at He, when the external field over- 
comes the exchange field. For applied field parallel to 
the hard axis the saturation occurs at H' E [14|. There 
is a spin-flop transition for field along the easy axis at 

H = h 'ea Il3 ■ 

The published susceptibility data, Fig. 4 of Ref. 

can be used to find all but one of the parameters of the 
model. The first peak in the dM/dH curve at 1.8 T has 
been already identified as the spin-flop field H' EA [5( . At 
around 8 T the dM/dH curve has a shoulder. We identify 
this as the onset of the saturation, and the corresponding 
field is He- Finally, the saturation is complete around 
13T, corresponding H' E . From these values we get A = 
8 T, a = 5 T and c = -0.42 T. 

The parameter b determines the ESR frequency 0J2 at 
zero field. We selected b = 6 T to match the zero field 
gap at 6.5 cm -1 (Mode (1) in Figures 01 and [5] This is 
quite close to the dispersionless magnon mode seen in 
neutron scattering (0.7 meV = 5.6 cm" 1 ) @- 

The principal frequencies are listed in Appendix B, and 



the results are shown in Fig. Considering that many 
of the parameters of the model were taken from other 
measurements, the agreement with the experiment is ex- 
cellent. The calculated zero field resonance uj\ ard is very 
close to the measured one at 53 GHz 5]. In fact, the low 
field behavior of our model (where the slope dw/dH is 
zero for H — > 0) fits the lower branch in Fig. 8 of Ref. 
|3 better than the uniaxial easy axis model used in that 
work (where the dio/dH = —7). Mode (1) is identified 
with u>2- As the field is increased, this mode is approx- 
imately independent of the field for any field direction. 
Therefore the powder signal remains narrow, and it is 
clearly visible in the fixed-field far-IR spectroscopy scans. 

Field-independent modes, like the io\ ard and w™* 
modes here, are very hard to detect in the field-sweep 
scan commonly used with the fixed-frequency methods. 
However, the spin-flop transition that occurs in the part 
of the sample where the field is close to the easy direc- 
tion causes a jump in the frequency of uf^ sy . This jump 
creates a measurable signal in the field-sweep scan |l3j . 
In our case the 225 GHz fixed frequency ESR measure- 
ment was just in the right frequency range to catch this 
feature. 

Above 5 T the powder average signal originating from 
L02 broadens and merges with other absorption. It can 
not be seen in the far-IR measurement, and only the 
(more sensitive) fixed frequency study picks it up as 
Mode (C). The upper extremum of all possible frequen- 
cies, , yields a strong signal. Although the g- factor 
was not adjusted (g — 2), the coefficient \Jl + H ci /He 
seen in Eq. ^] brings uJi right into the middle of the 
experimental values of the main branch (M). 

With the set of parameters A=8 T, a=5 T, b=6 T 
and c= -0.42 T, we can deduce the coupling energies in 
the Hamiltonian described by Eq. For the antifer- 
romagnetic coupling we get J AF — A/3— 3.58 K. The 
anisotropy of the ferromagnetic interaction is Jp — J F = 
(a + 6)/3=4.93 K, comparable to the antiferromagnetic 
coupling. The anisotropy of the antiferromagnetic cou- 
pling is weak J AF — J AF = (a — &)/3=0.45 K. The ferro- 
magnetic coupling, deduced from the Curie- Weiss tem- 
perature, J|,=-15 K, is the dominant interaction. This 
hierarchy of interactions was also deduced from inelas- 
tic neutron scattering by Lewis et a/.Q using a similar 
model. 

The weak anisotropy within the easy plane can be 
assigned to the ferromagnetic or antiferromagnetic cou- 
plings or a combination of the two. The parameter c rep- 
resents the sum of the two anisotropics: 2c/3 = (J F — 

Jf) ~ ~ Jf) + {Jaf ~~ ^af) ~ (^af ~ Jaf) = — 0.37 K. 
However, the measurements would be compatible with 
other combinations of these parameters (as long as the 
sum is fixed). For example all anisotropy may be in the 
ferromagnetic coupling {Jp — Jp) — (Jp — Jp) = —0.37 K 
and no anisotropy in the antifcrromagnetic coupling, 

(J%f ~ Jaf) ~ (Jaf ~ Jaf) = 0- 

We tested the model with various other choice of pa- 
rameters, partially disregarding the three characteris- 



tic fields obtained from the susceptibility measi 
Decreasing a to zero fails to describe the flatness 
(1). Negative values of a do not give good fits. '. 
the agreement with the data is still reasonabL 
special case the uniaxial anisotropy of the antif 
netic interplane coupling, J Z AF — J\ F , is exactb 

In spite of the general agreement between t. 
and the experiment, there are two features the 
unexplained: the splitting of mode (C) at low 
ture (Figure n insert) and the existence of moc 

The splitting if mode (C) is consistent with 01 
This mode is detected when the held is perpenc 
the hard axis. Note that the calculations for 
done with c — 0, but in the presence of a fin 
saturation field may be slightly different for stc 
along the easy and intermediate axes. The mag 
parameter |c| = 0.42 T is in the right range to explain 
the observed splitting. 

The existence of mode (2) is entirely beyond the simple 
model presented here. The two-sublattice Hamiltonian 
yields two modes at zero field, and we have those two 
modes identified as mode (A) at 1.75 cm^ 1 (53 GHz) 
and mode (1) at 6.5 cm -1 . No room is left for mode (2) 
in this model. 



V. CONCLUSION 

We studied the field dependence of the ESR in NaNiC>2 
and we interpreted the results in terms of a model involv- 
ing anisotropic exchange interactions. The main interac- 
tion is ferromagnetic in the triangular NiO plane. It is an 
order of magnitude larger than the interplane antiferro- 
magnetic coupling. The anisotropy in the ferromagnetic 
coupling is characterized by an uniaxial " easy plane" pa- 
rameter of a + b = 11 T, and a much weaker easy axis 
parameter of c = —0.42 T. 

One of the main results of our work is that the easy 
direction is only slightly better than the other directions 
perpendicular to the hard axis. This is in general agree- 
ment with the inelastic neutron scattering result, where 
the magnon spectrum was modeled with an easy plane 
spin Hamiltonian ||. For the easy plane system one of 
the magnon modes has zero energy (in zero field and at 
zero wavenumber). The anisotropy within the easy plane 
pushes this magnon mode to a finite energy (1.8 cm -1 = 
0.2 meV). This second magnon gap is small, and there- 
fore it is not visible in the neutron scattering, but it is 
clearly seen in the ESR studies. 

The principal axes of the anisotropy tensor were con- 
strained by symmetry considerations (see appendix A), 
but cannot be determined uniquely. In order to reconcile 
our results with the magnetization measurements on a 
twinned crystal and with the neutron diffraction re- 
sults 0, where the spin orientation was measured, the 
easy direction (the x axis of the reference frame) was se- 
lected so that it makes a 100° angle to the ab plane (see 
Fig. EI)- The hard axis was selected within the triangu- 




FIG. 6: The orientation of the hard axis z and the easy axis x 
relative to the triangular Ni layers. The ab plane is highlighted 
in the top layer. The z axis is parallel to the b direction in 
the crystal; the angle between the ab plane and the y axis is 
p « 10°. The spin directions are indicated by red arrows. 



lar nickel plane, along the b direction of the crystal (Fig. 
EJ. Alternatively, it may point in the y direction of the 
reference frame. The hard axis is probably related to 
spin-orbit effects on the ferromagnetic exchange. Indeed, 
the ratio \(J F — Jp/J F \ = 0.33 is comparable to the spin- 
orbit effect on the g tensor: (g 1 - — ,g")/.g" = 0.10. The 
choice of the easy axis involves an energy scale one order 
of magnitude smaller. 

The existence of a a third zero-field mode (mode 
2) may point beyond the model used here: any two- 
sublattice model, that assumes that the magnitude of 
the sublattice magnetisations are fixed, yields at most 
two resonance modes in zero field. One can then think 
of quantum effects. They may cause some modes to 
shift from zero to finite frequencies |l6j |. but do not 
change the number of modes. The comparison with 
other triangular magnets having antifcrromagnetic cou- 
pling within the layers is intriguing: (including CsNiBr3, 
RbNiBr 3 , CsMnI 3 , CsMnBr 3 CsCuCl 3 , RbCuCl 3 and 
RbFeQvloO^)- There, typically three or six ESR modes 
0,0,0 are observed. This is related to the umbrella- 
like spin order with, for a given triangle, each spin point- 
ing at 120 ° from each other. Three sublattices are def- 
initely present. This is not supposed to be the case in 
NaNi02 where there is strong support that the coupling 
within the triangular layers is ferromagnetic (the Curie- 
Weiss temperature is positive and the neutron diffraction 
patterns correspond to a A type antiferromagnet with fcr- 
romagnetically aligned spins in the layers [7J). However 
the sign of the interactions in the triangular planes has 
been the subject of intense debate in the isomorphic com- 
pound LiNi02 2]. Even for NaNi02 a recent calculation 
by Vernay et al |3j proposes that one out of three inter- 
actions bond in the triangles are antiferromagnetic lead- 
ing to antiferromagnetic chains, although this magnetic 
structure has not been observed experimentally. More 
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over, in the intermediate compound Lio.3Nao.7Ni02 @, 
which, on the orbital sector, behaves as LiNiC>2 and, on 
the spin sector, presents a long range antiferromagnetic 
order as in NaNiC>2, we also observe the occurrence of 
three modes at zero field (around 0, 5.2 and 6.5 enr 1 ); 
there the magnetic structure has not been solved yet. It 
is tempting to assign these similar behaviors to the pres- 
ence of additional antiferromagnetic interactions. Fur- 
ther studies, preferably on single crystals, may solve this 
remaining puzzle in the NaNi02/LiNi02 layered trian- 
gular compounds. 
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VII. APPENDIX A: DERIVATION OF THE 
MEAN FIELD HAMILTONIAN AND 
SYMMETRIES 

In the molecular field approximation the spins are as- 
signed to two sublattices: Spins in the even layers are 
labeled by \i and spins in the odd layers are labeled by v. 
Si and Sj in the Hamiltonian are replaced with their av- 
erage values, Si = 2/N and S 2 = 2/N S^, where 
N is the total number of spins. 

First, let us look at the term ^2 ab Si J a bSj in the 
Hamiltonian Eq. 0] All of the spins in this sum belong 
to the same sublattice, so the terms to add up will have 
the form of Si J ;,Si and S2 J a bS 2 . For any given spin 
within a ferromagnetic layer, there are six neighbors, ap- 
proximately positioned on the corners of a hexagon. Even 

though the exchange J a \, with any one of this neighbors 
may be anisotropic, in the sum of the six terms most of 



the anisotropy within the plane will be cancelled. In fact, 
as long as the exchange couplings can be represented by 
tensors, in an undistorted lattice the cancellation will be 
exact. The corresponding part of the Hamiltonian would 
be: EafcS/j afc S, = (N/2) (l/2)6(Si JfSi + S 2 JfS 2 ), 

where the effective coupling Jp is isotropic within the 
plane, the extra factor of 1/2 compensates for the double 
counting of each bond and the factor 6 reminds us that 
this exchange is an average over six bonds. Similar ar- 
guments can be used for the antiferromagnetic interlayer 
coupling to derive J af from a sum over J c 's. The 
proper reference frame for the new coupling tensors has 
the x direction perpendicular to the layers. The (equiv- 
alent) y and z directions are within the layers. 

The anisotropy of the effective couplings J p and 
J af is related to the distortion of the lattice. While 
above 480 K all nearest neighbor Ni-Ni bond lengths are 
equal 0, below the transition one side of the triangu- 
lar lattice (along the b direction) is compressed, and two 
other sides are elongated. The a — c plane remains a mir- 
ror plane, but the direction perpendicular to the layers 
is not a high symmetry direction anymore. (See Fig. 0) 
Accordingly, the x direction will tilt relative to the direc- 
tion perpendicular to the layers, as it acquires a non-zero 
component in the crystallographic a direction. At the 
same time the anisotropy axes within the plane become 
well defined. We will take the z axis parallel to the b crys- 
tallographic direction. The y axis will point slightly out 
of the plane, being perpendicular x and z. Notice that, 
except for permutations of the x,y,z labels, this is the 
only selection satisfying the broken crystal symmetries. 
The magnitude of the tilt of the x axis, and the amount of 
anisotropy, remain a free parameters. Similar arguments 
apply to the nearest neighbor antiferromagnetic coupling 
between the layers. 

Spin resonance on a powder sample can be used to eval- 
uate the anisotropy parameters of the resonance line. As 
we will see later, our system is described by a model with 
a dominant " easy plane" anisotropy, and a much weaker 
anisotropy within the plane. However, the orientation of 
the anisotropy axes relative to the crystal cannot be can 
not be determined from a powder measurement. Neu- 
tron diffraction 7] and single crystal magnetization data 
indicate that the easy axis makes an angle of approxi- 
mately 100° to the NiO planes. We select this to coincide 
with the easy direction (the x axis of the reference frame) . 
The hard axis may be either along y or z; both choices 
agree with existing symmetries. We tentatively pick the 
z direction for the hard axis. 

The mean field Hamiltonian becomes: 
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Umf - — [J F (Sf x + S$ x ) + Jp(Sf y + Sly) + J z F {Sl z + S 2 Z )] + 

+3iV (J% F S lx S 2x + J y AF S ly S 2 y + J% F S u Sl) + y<?Ms(Si + S 2 )H (6) 

The dimensionless sublattice magnetizations are denned as mi = ^Si/S and 1112 = ^S 2 /S. When thermal and 
quantum fluctuations are neglected, the length of these vectors is fixed | mi | = | m 2 \ = 1/2. Accordingly, we can express 
one of the sublattice magnetization components with the other two: m\ z — \/A~m\ x —m\ y &ndm 2 z = 1/4— m\ x —m\ y . 
The Hamiltonian is: 



Hmf = 12NS 2 [{J% - J F ){m\ x + m\ x ) + (J* - J F )(mj y + m\ y )\ + 
+6NS 2 (J AF m lx m 2x + J v AF rn\ y m 2y + JAF m izm 2z ) + NSg(i B (mi + m 2 )H (7) 

and we dropped a constant term proportional to J p. In the absence of fluctuations the magnitude of the ferromagnetic 
coupling becomes irrelevant, only the anisotropy matters. Interestingly, the anisotropy term obtained here for the 
ferromagnetic coupling cannot be distinguished (within the molecular field approximation) from a microscopic " single 
ion" spin anisotropy. However, the physical origins are different: For spin 1/2 the microscopic single ion anisotropy 
can be shown to be exactly zero, whereas the same arguments do not apply here, where the anisotropy is due to the 
ferromagnetic exchange coupling. 

The easy plane anisotropy around the z axis corresponds to (Jp — Jp) « (Jp — Jp) (and similarly for the antifer- 
romagnetic coupling). We can re-arrange the terms to better reflect this fact. For mi we get 

(J F - J F )m\ x + (J F - J F )m% = 
= 1/2[(J| - J F ) + (J| - J F )](m 2 lx + m\ y ) + \/2[(J F - J F ) - (J F - J F )](m\ x - m\ y ) = 
= 1/2[(J F - J F ) + (J F - J F )](1 - ml) + 1/2[(J F - J F ) - (jy - J F )](2 m \ x - 1 + m 2 u ) = 

- (Jp - J F )m\ z + [(Jp - J F ) - (J F - J F )\m\ x + const. (8) 

This way the coefficient of the m\ z term represents the uniaxial anisotropy; the (smaller) coefficient of the m\ x term 
describes the remaining anisotropy around the z axis. Similar rearrangement can be made for 7712 and for all of the 
antifcrromagnetic terms. 

The total magnetization and the antiferromagnetic order parameter vectors are introduced as m = (mi + 1112) 
and 1 = (mi — m2). Notice that |m| = 1 corresponds to the saturation magnetization, and |1| = 1 describes perfect 
antiferromagnetic order. In terms of these new parameters the Hamiltonian can be re-written as: 



Hmf = 3NS 2 [J z AF (m 2 - l 2 ) + (J AF - J v AF )(m 2 z - l 2 z ) + (J% F - J* AF ) - (J\ F - J z AF )(m 2 x - l 2 x ) 

+ (J F ~ J V F )(ml + ll) + +[(J x f ~ Jf) ~ (J V f ~ J¥)« + ID] + NgfismH (9) 

Here we re-grouped the terms in the antiferromagnetic coupling, emphasizing the anisotropy, without introducing new 
constrains. Neglecting fluctuations leads to ml = and m 2 + 1 2 = 1. Eliminating the l 2 term finally yields: 

Hmf/N = ^m 2 + %n 2 z + b -l 2 z + C -m 2 x + U 2 X - mh (10) 

where A = 12S 2 J AF , h = g^sST!, and the anisotropy constants are expressed as 

a = 6S 2 (J F — Jp + J Z AF — J AF ), b = 6S 2 (Jp — Jp — J AF + J AF ), 

c = 6S 2 [(J F - J F ) - (J F J F ) + (J\ F - J Z AF ) - (J\ F - J Z AF )\ 

d = 6S 2 [(J F - J Z F ) - (J F - J Z F ) - (J\ F - J Z AF ) + (J\p - J Z AF % (11) 

Notice that Eq. ^|is applicable in the most general case, but the coupling constants reflect an expected hierarchy. 
A simple isotropic antiferromagnet is obtained for a = b = c = d = 0. For uniaxial anisotropy c = d = 0; the 
coefficients a and b describe anisotropies in the total magnetization and antiferromagnetic order, respectively. In our 
case coefficients c and d will be small. , 



VIII. APPENDIX B: SPIN RESONANCE gleet the small anisotropy within the easy plane, and take 

FREQUENCIES. 



For any given field direction there are two modes, de- 
noted by uji and lo 2 - At high external fields one can ne- 
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c = d = 0. The frequencies were determined by Turov 
|14|. There are two principal frequencies for each modes, 
corresponding to external field applied parallel and per- 
pendicular to the hard axis. The results are expressed in 
terms of effective fields: 

H E = A/M 
H' E = (A + a)/M 
H E = (A-b)/M 
H a = a/M 
H b = b/M 



He a = v HsHb 



H \\ - V H E H b 

E 



where Mq is the saturation magnetization. 

For field perpendicular to the hard axis one gets: 



(12) 



u>i = jHy/l + H a j H E if H <H E 
= jy/ H(H + H a ) if H > H E 

and 

ui= 1 H EA Jl-H*/H% 



if H < H E 



iJ(H-H%)(H-He) if H>H£13) 



where 7 = g/is/Ti in terms of the Bohr magneton and 
the g-factor. He is the saturation field; for H > He all 
spins line up with the external field. 

For the field direction parallel to the hard axis the 
saturation field is H' E > He, since both the exchange 
coupling and the anisotropy work against the external 
field. The frequencies are: 



,11 



if 



H<H' e 
if 



H>H' E 



and 



4 = iHea^I + H 2 /H\ if H < H' E 
4 = j(H - H a ) if H > H' E 



(14) 



The modes are shown in the high field part of Fig. [21 
The u>2 mode reaches zero at the lower saturation field 
He- Below the saturation field the lo^ mode looks like 
a "free" spin resonance, except for the apparent g-factor 
is increased by a factor of \Jl + H a /HE- Well above the 
saturation field the line is shifted to higher frequencies 
by the amount of jH a /2, but its slope is still 7. For 

fields parallel to the hard axis the uj\ mode is zero at 
the upper saturation field. With the particular choice 

of parameters indicated in the Figure. The oj\ mode is 
nearly independent of the field up to H' E . 



At low fields the anisotropy within the easy plane be- 
comes important. A finite value of the parameter c (or 
d) has two important consequences. First, in zero field, 
the 0J1 mode (the "Goldstone mode") will shift to finite 
frequency. Second, the field dependence of the uj^ mode 
will be different for external fields applied in different 
directions within the easy plane. Instead of the two prin- 
cipal directions relative to the hard axis z (parallel and 
perpendicular), we will have to deal with three axes. We 
will identify these directions as the hard axis (parallel to 
z), the intermediate axis and the easy axis (formerly the 
two equivalent directions perpendicular to z). 

For an approximate treatment of the c ^ case we 
will turn to an early paper by Nagamiya |13| . assuming 
// ///;. 

In finite fields along the hard axis the u>i mode is ap- 
proximately independent of the field (this mode corre- 
sponds to the rotation of the spins in the easy plane, and 
it is excited with an oscillating magnetic field parallel to 
the static field). For fields in the intermediate direction 
the frequency of the mode increases; for the field applied 
in the easy direction there will be a spin-flop transition 



at H = m 



EA 



'Ac I 'M - 



Jiard 



= j'H 



EA 



EA 



\ - + H 2 - v//? 2 /4 + 2aH 2 if H<H' t 



L EA 



EA 



(15) 



Here we used a — H E a 



TT/2 

H EA 



T 2 

l EA 



TT/2 

H EA- 



(3 = Hi 

Nagamiya's calculation does not account for the increase 
of the apparent g-factor seen in Turov's low field result. 
Therefore we are using a modified 7' = 7-^/1 + H a j He 
here. 

The anisotropy within the easy plane has no influence 
on the ujlj mode (for fields parallel to the hard axis) . For 
fields applied in the intermediate direction, the u>2 mode 
is approximately independent of the field (except for a 
slight decrease in frequency, a precursor of the eventual 
saturation at He)- When the field is in the easy direc- 
tion, the spin-flop transition results in a jump of the ESR 
frequency. 



hard 



= lH, 
1 H J 



EA 



EA 



= 7-1/77 +H 2 + V/3 2 /4 + 2aH 2 if H<H' t 



easy 
1 y — 



1X /H 2 EA - H' 2 A 



if H > H' f 



EA 



EA 



(16) 



To ensure a seamless match to the high field calcula- 
tion, in Fig. 0]we multiplied these frequencies with a sat- 
uration term, yfl — H 2 / H E . For the hard axis direction 
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we used the original w\ values. Neither this procedure, trary fields and parameters, but they are qualitatively 
nor the modified 7' introduced earlier, can completely correct, 
substitute for the full solution of the problem for arbi- 
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